On the Spectrum of a Quantum Dot 
with Impurity in the Lobachevsky Plane 

P. Stovicek and M. Tusek 

Abstract. A model of a quantum dot with impurity in the Lobachevsky plane 
is considered. Relying on explicit formulae for the Green function and the 
Krein Q— function which have been derived in a previous work we focus on 
the numerical analysis of the spectrum. The analysis is complicated by the fact 
that the basic formulae are expressed in terms of spheroidal functions with 
general characteristic exponents. The effect of the curvature on eigenvalues 
and eigenfunctions is investigated. Moreover, there is given an asymptotic 
expansion of eigenvalues as the curvature radius tends to infinity (the flat 
case limit). 
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1. Introduction 

The influence of the hyperbolic geometry on the properties of quantum mechanical 
systems is a subject of continual theoretical interest for at least two decades. 
Numerous models have been studied so far, let us mention just few of them [1, 2, 
3, 4]. Naturally, the quantum harmonic oscillator is one of the analyzed examples 
[5, 6]. It should be stressed, however, that the choice of an appropriate potential 
on the hyperbolic plane is ambiguous in this case, and several possibilities have 
been proposed in the literature. In [7], we have modeled a quantum dot in the 
Lobachevsky plane by an unbounded potential which can be interpreted, too, as 
a harmonic oscillator potential for this nontrivial geometry. The studied examples 
also comprise point interactions [8] which are frequently used to model impurities. 

A Hamiltonian describing a quantum dot with impurity has been introduced 
in [7]. The main result of this paper is derivation of explicit formulae for the 
Green function and the Krein Q— function. The formulae are expressed in terms of 
spheroidal functions which are used rather rarely in the framework of mathematical 
physics. Further analysis is complicated by the complexity of spheroidal functions. 
In particular, the Green function depends on the characteristic exponent of the 
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spheroidal functions in question rather than directly on the spectral parameter. 
In fact, it seems to be possible to obtain a more detailed information on eigenvalues 
and eigenfunctions only by means of numerical methods. The particular case, when 
the Hamiltonian is restricted to the eigenspace of the angular momentum with 
eigenvalue 0, is worked out in [9]. In the current contribution we aim to extend the 
numerical analysis to the general case and to complete it with additional details. 

The Hamiltonian describing a quantum dot with impurity in the Lobachevsky 
plane, as introduced in [7], is a selfadjoint extension of the following symmetric 
operator: 

H = - + I C °0 I + h Sinh " (o) W + 4a 5 ) + I aW Sinh2 (!) ' 

Dom(iJ) = C^°((0,oo) x S 1 ) C £ 2 ((0,oo) x S 1 , a smh(g/ a)dg d<t>) , 

where (g, <j>) are the geodesic polar coordinates on the Lobachevsky plane and a 
stands for the so called curvature radius which is related to the scalar curvature 
by the formula R — —2/ a 2 . The deficiency indices of H are known to be (1,1) 
and we denote each selfadjoint extension by H(x) where the real parameter \ 
appears in the boundary conditions for the domain of definition: f(g,4>) belongs 
to Dom(H(x)) if there exist /o, fi 6 C so that /i : /o = X : 1 an d 

f(g, <j>) = ~fo log(g) + fi + o(l) as g - 0+ 

(the case x = oo means that / = and /i is arbitrary), see [7] for details. H(oo) is 
nothing but the Friedrichs extension of H . The Hamiltonian H(oo) is interpreted 
as corresponding to the unperturbed case and describing a quantum dot with no 
impurity. 

After the substitution £ = cosh(g/a) and the scaling H = a~ 2 H, we make use 
of the rotational symmetry (which amounts to a Fourier transform in the variable 
</>) to decompose H into a direct sum as follows 

oo 

H= H m , 

m— — oo 

Dom(H m ) = Cg°(l, oo) C L 2 ((l, oo), d^). 

Let us denote by H m , m e Z, the restriction of H(oo) to the eigenspace of 
the angular momentum with eigenvalue m. This means that H m is a self-adjoint 
extension of a~ 2 H m . It is known (Proposition 2.1 in [7]) that H m is essentially 
selfadjoint for ra/0. Thus, in this case, H m is the closure of a~ 2 H m . Concerning 
the case m = 0, H is the Friedrichs extension of a~ 2 Ho- For quite general reasons, 
the spectrum of H m , for any to, is semibounded below, discrete and simple [10]. 
We denote the eigenvalues of H m in ascending order by E n , m (a 2 ), n e No- 
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The spectrum of the total Hamiltonian H{\), x ^ 00, consists of two parts 
(in a full analogy with the Euclidean case [11]): 

1. The first part is formed by those eigenvalues of H(x) which belong, at the 
same time, to the spectrum of if (00). More precisely, this part is exactly the 
union of eigenvalues of H m for m running over Z \ {0}. Their multiplicities 
are discussed below in Section 5. 

2. The second part is formed by solutions to the equation 

Q H (z) = x (i.i) 

with respect to the variable z where Q H stands for the Krein Q-function of 
iJ(oo). Let us denote the solutions in ascending order by e„(a 2 ,x), n G No- 
These eigenvalues are sometimes called the point levels and their multiplicities 
are at least one. In more detail, e n (a 2 ,x) is a simple eigenvalue of H{\) if it 
does not lie in the spectrum of H(oo), and this happens if and only if e„(a 2 , x) 
does not coincide with any eigenvalue Ei^ m (a 2 ) for I £ Nq and meZ,ra/0. 

Remark. The lowest point level, eo(a 2 , x), lies below the lowest eigenvalue of H(oo) 
which is Eofi(a 2 ), and the point levels with higher indices satisfy the inequalities 
£„_i,o(a 2 ) < e«(a 2 ,x) < E nfl (a 2 ), n= 1,2,3,... . 



2. Spectrum of the unperturbed Hamiltonian H(oc) 

Our goal is to find the eigenvalues of the rath partial Hamiltonian H m , i.e., to find 
square integrable solutions of the equation 

or, equivalently, 

H m ^,(0 = a 2 z^(i). 

This equation coincides with the equation of the spheroidal functions (A.l) pro- 
vided we set fi = \m\, 6 = — a 4 w 2 /16, and the characteristic exponent v is chosen 
so that 

" V 16 J 4 
The only solution (up to a multiplicative constant) that is square integrable near 
infinity is sj, m|(3) (£, -a 4 w 2 /16). 

Proposition A. 3 describes the asymptotic expansion of this function at £ = 1 
for m G N. It follows that the condition on the square integrability is equivalent 
to the equality 

^+1/2)^ _,(-^ +K ?(-^ =0. (2.1) 

Furthermore, in [7] we have derived that 

(£, 6) = a log(£ -I) +13 + 0((£ - 1) log(£ - 1)) as £ -> 1+, 
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where 



a = 



itan(^7r) e 



-i(2v+l/2)ir 



Taking into account that the Friedrichs extension has continuous eigenfunctions 
we conclude that equation (2.1) guarantees square integrability in the case m = 0, 
too. 



FIGURE 1. Eigenvalues of the partial Hamiltonian Hi 




As far as we see it, equation (2.1) can be solved only by means of numerical 
methods. For this purpose we made use of the computer algebra system Mathe- 
matica 6.0. For the numerical computations we set w = 1. The particular case 
m = has been examined in [9]. It turns out that an analogous procedure can 
be also applied for nonzero values of the angular momentum. As an illustration, 
Figure 1 depicts several first eigenvalues of the Hamiltonian H\ as functions of 
the curvature radius a. The dashed asymptotic lines correspond to the flat limit 
(a — > oo). 

Denote the nth normalized eigenfunction of the mth partial Hamiltonian H m 
by ip n,m{£)- Obviously, the eigenfunctions for the values of the angular momen- 
tum m and — m are the same and are proportional to Si m '^(£, — a 4 o> 2 /16), with 
v satisfying equation (2.1). Let us return to the original radial variable g and, 
moreover, regard H m as an operator acting on L 2 (R + ,dg). This amounts to an 
obvious isometry 

L 2 (M+, a' 1 smh(g/a)dg) -► L 2 (M+, dg) : f(g) ^ a' 1 ' 2 smh 1/2 (g/a)f(g). 

The corresponding normalized eigenfunction oiH m , with an eigenvalue a 2 z, equals 

(\ \ 1/<2 
ipn,m{Q) = ( - sinh (J ) J Vi„ ;m (cosh . (2.2) 
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At the same time, relation (2.2) gives the normalized eigenfunction of H m (consi- 
dered on L 2 (R + ,dg)) with the eigenvalue z. The same Hilbert space may be used 
also in the limit Euclidean case (a = oo). The eigenfunctions $ n , TO in the flat case 
are well known and satisfy 



The fact that we stick to the same Hilbert space in all cases facilitates the 
comparison of eigenfunctions for various values of the curvature radius a. We 
present plots of several first eigenfunctions of Hi (Figures 2, 3, 4) for the values 
of the curvature radius a = 1 (the solid line), 10 (the dashed line), and oo (the 
dotted line). Again, see [9] for analogous plots in the case of the Hamiltonian H . 
Note that, in general, the smaller is the curvature radius a the more localized is 
the particle in the region near the origin. 



3. The point levels 

As has been stated, the point levels are solutions to equation (1.1) with respect to 
the spectral parameter z. Since, in general, Q(z) = Q(z) the function Q(z) takes 
real values on the real axis. Let H(oo) = a 2 H(oo) be the Friedrichs extension of 
H. An explicit formula for the Krein Q-function Q H (z) of H(oo) has been derived 



Figure 2 . The first eigenfunction of the partial Hamiltonian Hi 

i/<o,i(p)fora 2 =l,10 ,00 
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in [7]: 

^) = -^(- log(2 )- 2 , (1) + 2 ^(-^) S o(-^)) 

, 1 ( p «r(3„ + 3/2) _ , V 1 , (2^) 

2a 2 tan(^7r) ^ #0(_<il|i) J + 4?ra 2 

where 1/ is chosen so that 



-z . 



16 J 4 
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The symbol K^(0) stands for the so called spheroidal joining factor, 

oo 

*8 V (0):= ]T (-l) r a° Vtr (6)*{u + l + 2r), 

r— — oo 

where the coefficients a° r (0), r G Z, come from the expansion of spheroidal func- 
tions in terms of Bessel functions (for details see [7, the Appendix])), and s®(6) is 
defined by formula (A. 4). One can obtain the Krein Q-function of H(oo) simply 
by scaling Q H (z) = a 2 Q H {a 2 z). 




Since we know the explicit expression for the Krein Q-function as a function 
of the characteristic exponent v rather than of the spectral parameter z itself 
it is of importance to know for which values of v the spectral parameter z is 
real. Propositions A.l and A. 2 give the answer. For i/£l and for v of the form 
v = —1/2 + it where t is real, the spheroidal eigenvalue A™(— a 4 w 2 /16) is real, 
and so the same is true for z. Moreover, these values of v reproduce the whole 
real z axis. With this knowledge, one can plot the Krein Q-function Q H = Q H (z) 
for an arbitrary value of the curvature radius a. Note that for a — oo, the Krein 
Q-function is well known as a function of the spectral parameter z [12] and equals 
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(setting uj = 1, ^ is the logarithmic derivative of the gamma function) 
Q(z) = ^ (-* (^) + log(2) + 2tt(l)) . 

Next, in Figure 5, we present plots of the Krcin Q-function for several dis- 
tinct values of the curvature radius a. Moreover, in Figure 6 one can compare the 
behavior of the Krein Q-function for a comparatively large value of the curvature 
radius (a 2 = 24) and for the Euclidean case (a — oo). 

Figure 6. Comparison of the Krein Q-functions for a 2 = 24 and 

a 2 = oo 

a 2 =24, a 2 =oo (dashed) 




z 



Again, equation (1.1) can be solved only numerically. Fixing the parameter 
X one may be interested in the behavior of the point levels as functions of the 
curvature radius a. See Figure 7 for the corresponding plots, with \ = 0, where 
the dashed asymptotic lines again correspond to the flat case limit (a = oo). Note 
that for the curvature radius a large enough, the lowest eigenvalue is negative 
provided x is chosen smaller than Q(0) ~ 0.1195. 

4. Asymptotic behavior for large values of a 

The mth partial Hamiltonian H m , if considered on L 2 (R + ,dg), acts like 

° 2 m2 - I 1 2 2 • , 2(Q\ d 2 . , 

H m = -^-o + 9 . ; 2 , as + 7^ smh - = : ~7T~2 + v m{a, en- 
ds 2 a 2 smli (§) 4 W dg 2 

For a fixed g ^ 0, one can easily derive that 

ym(a,g) = -^ + -^g 2 + ^ r + Y ^ + 0^- 4 j asa^oo. 
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Figure 7. Point levels for H(0) 



e ; (a 2 , 0), 1=0,1,2,3 




Recall that the mth partial Hamiltonian of the isotropic harmonic oscillator on 
the Euclidean plane, if considered on L 2 (R + ,dg), has the form 



dg< 



+ 



m 



- + r v - 



This suggests that it may be useful to view the Hamiltonian H m , for large values 
of the curvature radius a, as a perturbation of 



+ - 4m 2 + c 2 /) =: H% + -^U m (g). 



12a 2 v * ' ' m 12a 2 

The eigenvalues of the compared Hamiltonians have the same asymptotic expan- 
sions up to the order 1/a 2 as a — > oo. 

Let us denote the nth eigenvalue of the Hamiltonian by E^ m , n G No- 
It is well known that 

E% m = (2n+\m\ + l)u> 

and that the multiplicity of E^ m in the spectrum of H E equals 2n + \m\ + 1. The 
asymptotic behavior of E n ^ m (a 2 ) may be deduced from the standard perturbation 
theory and is given by the formula 



E n ,m{a 2 ) — E^ r 



1 {^n,m^^m^n,m} 

12a 2 ($„, m ,$„ .ml 



as a 



oo, 



(4.1) 



where $ n ,m denotes a (not necessarily normalized) eigenfunction of H E associated 
with the eigenvalue E^ m (see (2.3)). The scalar products occurring in formula 
(4.1) can be readily evaluated in L 2 (R + , dg) with the help of Proposition A. 4. The 
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Table 1. Comparison of numerical and asymptotic results for 
the eigenvalues, a 2 = 24 

£-o,o -£q,o i?2,o Ep y i E\ t \ E2S 

numerical 1.0265 3.162 5^42 2.060 4.259 6.58 

asymptotic 1.0268 3.169 5.46 2.058 4.258 6.59 

error (%) -0.03 -0.22 -0.74 0.10 0.02 -0.15 



resulting formula takes the form 

E n , m {a 2 ) = (2n + |ro| + 1) w + ^2n(n + \m\ + 1) + |m| + ^ + o(J^ (4.2) 

as a — ► 00. This asymptotic approximation of eigenvalues has been tested nume- 
rically for large values of the curvature radius a. The asymptotic eigenvalues for 
a 2 = 24 are compared with the precise numerical results in Table 1. It is of interest 
to note that the asymptotic coefficient in front of the a~ 2 term does not depend 
on the frequency u>. 



5. The multiplicities 

Since H- m = H m the eigenvalues E n _ m (a 2 ) of the total Hamiltonian H(oo) are at 
least twice degenerated if m ^ 0. From the asymptotic expansion (4.2) it follows, 
after some straightforward algebra, that no additional degeneracy occurs and thus 
theses eigenvalues are exactly twice degenerated at least for sufficiently large values 
of a. 

Applying the methods developed in [11] one may complete the analysis of the 
spectrum of the total Hamiltonian H(x) for \ 00 ■ Namely, the spectrum of H{x) 
contains eigenvalues E n ^ m (a 2 ), m > 0, with multiplicity 2 if Q H (E n ^ m (a 2 j) ^ x, 
and with multiplicity 3 if Q H (E n , m (a 2 )) = \- The rest of the spectrum of H(x) is 
formed by those solutions to equation (1.1) which do not belong to the spectrum of 
H(oo). The multiplicity of all these eigenvalues in the spectrum of H{\) equals 1. 



Appendix: Auxiliary results 

In this appendix we summarize several auxiliary results. Firstly, for our purposes 
we need the following observations concerning spheroidal functions. The spheroidal 
functions are solutions to the equation 

(! - - 2^ + [K{0) + 40(1 e) M 2 (l - a" 1 ] * = 0. (A.l) 

For the notation and properties of spheroidal functions see [13]. A detailed infor- 
mation on this subject can be found in [14], but be aware of somewhat different 
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notation. A very brief overview of spheroidal functions is also given in the Appen- 
dix of [7]. 

In the last named source, the following proposition has been proved in the 
particular case m = 0. But, as one can verify by a direct inspection, the proof 
applies to the general case m E Z as well. 

Proposition A.l. Let v,6 el, m eZ. Then A™(6>) E E. 

The following claim is also of interest. 

Proposition A. 2. Let v = —1/2 + it where t E E, and 9 E E, m E Z. Then 
A™(6>) G E. 

Proof. Let us recall that the coefficients a™ r (0) in the series expansion of spheroidal 
functions in terms of Bessel functions satisfy a three term recurrence relation (see 
[7, the Appendix]), 

^r(«Kr-l(«) + <r(«)<rW + ^Kr+M = WKrW- (A-2) 

One may view the set of equations (A. 2), with r E Z, as an eigenvalue equation 
for A™ (6*) that is an analytic function of 8. A particular solution is fixed by the 
condition A™(0) = v(y + 1). Consider the set of complex conjugated equations. 
Since pf^ = (3^ r {6) = l3 m „_ l r {0), and the similar is true for a™ r (0) and 7™ r (0), 
it holds true that 



n-i, r WC-iW+^-i, r ^W+7™-i, r W< r+ iW = ^(«)< r (0). 
Since for each v of the considered form, 



A™ _i(0) = (-1/ - !)(-!/) = 1/(1/ + 1) = + 1) = A-(0), 



one has \™ v _ x {6) = A™(0). Moreover, \ m v _ 1 {6) = A™(0) in general. We conclude 
that \™{0) E E. ' □ 

Another auxiliary result concerns the asymptotic expansion of the radial 
spheroidal function of the third kind. 

Proposition A.3. Let v £ {-1/2 + fc| k E Z} , m E N. Tftera 

^'(6*) na?(0)e-«»+W* [ K —^ + J K " D 

as £ -» 1 + • (A.3) 
Proof. By the definition of the radial spheroidal function of the third kind, 

^C (3) M) := ■ 7 , fe^M) +ie-*^ Sr (1) (£,0)) , 

ZCOS(l/7Tj V / 

and by the relation between the radial and the angular spheroidal functions, 
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one has 

Using the definition 

oo 

Qs?(Z,0)= E (-i) r <r(<W +2r (0 

r— — oo 

and due to the well known asymptotic expansion for the Legendre functions [13], 

QT(0 ~ (-l) m 2 ro / 2 - 1 r(m)(^ - l)- m ' 2 as £ - 1+ , 
one derives that 

where 

oo oo 

~ E (-!)XrW = E (-l) r a-„-i,-r(0) = (C-iW)" 1 - 

r— — oo r— — oo 

(A.4) 

Hence <3s™j,_i(£, #) ~ Qs™(£, 0) as £ — > 1+, and one immediately obtains (A. 3). 

□ 

Further some auxiliary computations follow that we need for evaluation of 
scalar products of eigenfunctions (see (4.1)). 

Proposition A.4. Let \Fi{a,h,t) stand for the Kummer confluent hypergeometric 
function, andn,m,l G No- Then 

/■OO 

/ i m+, e-*iFi(-n,l + m,t) 2 dt 
Jo 

2 ^ +k (n\ (fc + Q! /fc + m + A ( A ' 5 ) 

fc=max{0,n-0 V 7 V ; V 7 

Proof. By definition, 



n / \ j-k n / 

■'-(-•i+».')-g(^ H -™ig(-')'(; 



kj (m + k)\ ' 



Let us denote the LHS of (A. 5) by /. Then the integral representation of the 
gamma function implies 

k( n \ ( n \ (j + k + m + l)\ 

j,k=0 



I = (ml? E (-!) j+fe U w n • ( A - 6 ) 

^ v Vj/ W (m + j)!(m + fc)! 
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Partial summation in (A. 6) can be carried out, 




xY) 



x=l 



(A.7) 



Expression (A.7) vanishes for k < n — I and equals 



{-iY(k + i)\( k + m + l \ 

\ n + m J 



for k > n — I. The proposition follows immediately. 



□ 



Corollary A. 5. In the case I = 0, (A.5J takes a particularly simple form: 
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